Rashba spin torque in ferromagnetic metal films 
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In a two-dimensional ferromagnetic metal film lacking inversion symmetry, the itinerant electrons 
mediate the interaction between the Rashba spin-orbit interaction and the ferromagnetic order 
parameter, leading to a Rashba spin torque exerted on the magnetization. Using Keldysh technique, 
in the presence of both magnetism and spin-orbit coupling, we derive a spin diffusion equation that 
provides a coherent description to the diffusive spin dynamics. The characteristics of the spin torque 
and its implication on magnetization dynamics are discussed for a wide range of relative strength 
between Rashba spin-orbit coupling and the ferromagnetic exchange. 
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I. INTRODUCTION 

By transferring angular momentum between electronic 
spin and orbital degrees of freedom, spin-orbit coupling 
fills the need for electrical manipulation of spins. Out- 
standing examples are the electrically generated bulk 
spin polarization-^ and the well-known spin-Hall effect 
(SHE) 3-5 in a two dimensional electron gas where the 
spin-orbit interaction, particularly of the Rashba-type,^ 
plays the leading role. Rashba spin-orbit interaction not 
only introduces an effective field perpendicular to the lin- 
ear momentum but also provides the backbone to the 
spin relaxation through the so-called D'yakonov-Perel 
mechanism^ which is dominant in a two-dimensional 
system. Besides its prominent role in semiconductors, 
Rashba spin-orbit coupling is believed to exist at fer- 
romagnetic/heavy metal and ferromagnetic/metal-oxide 
interfaces, where the inversion symmetry breaking offers 
a potential gradient empowering the spin-orbit coupling. 

Meanwhile, magnetism continuously stimulates the in- 
dustrial and academic appetite. In the pursuit of fast 
magnetization switching, Slonczewski-Berger spin trans- 
fer torque^ employs a polarized spin current instead of 
a cumbersome magnetic field. This celebrated scheme 
demands noncollinear magnetic textures in forms of, for 
example, spin valves or domain wall structures: 9 

In the presence of inversion symmetry breaking (such 
as asymmetric interfaces), a ferromagnetic metal layer as- 
sembles both magnetism and spin-orbit coupling, offering 
an alternative switching mechanism >i2r— spin-orbit cou- 
pling transfers the orbital angular momentum carried by 
an electric current to the electronic spin, thus creating 
an effective magnetic field (Rashba field). As long as the 
effective field is misaligned with the magnetization direc- 
tion, the so-called Rashba torque emerges to excite the 
magnetization. 

Current-driven magnetization dynamics by spin-orbit 
torque has been demonstrated by several experiments 
in metal oxide based systems*^— In fact, the Rashba 
torque can be categorized into a broader family of 
spin-orbit interaction-induced torques that has been 
observed in diluted magnetic semiconductors *£r— Re- 



cently, Miron et al.^ has demonstrated the current- 
induced magnetization switching using a single ferro- 
magnet in Pt/Co/AlOa; trilayers, which further consol- 
idates the feasibility of Rashba torque. The same type 
of torque is predicted to improve current-driven domain 
wall motion , 11 ! 20 which is also supported by experimen- 
tal observations.— At this stage, we are mindful of an 
alternative explanation based on the spin-Hall effect oc- 
curring in the underlying heavy metal layer (such as Pt or 
Ta)i22 An in-depth discussion on the distinction between 
the spin-Hall effect induced torque (called SHE torque 
thereafter) and the Rashba one is offered elsewhere^ 3 - 

In searching for a general form of the Rashba torque 
in an ultrathin ferromagnetic metal film, we found an 
expression that consists of two components J^an in-plane 
torque (aligned along the direction m x (y x m)) and 
an out-of-plane one (aligned to y x m), given y is the 
in-plane direction transverse to the current transport di- 
rection and m is the orientation of magnetization. In 
the case of a weak Rashba spin-orbit coupling, numerical 
solutions on a two-dimensional nanowire with one open 
transport direction has been carried out to appreciate the 
significance of diffusive motion on the spin torque. We 
found that the magnitude of the in-plane Rashba torque 
is enhanced by narrowing the magnetic wireji 2 - 

In this article, we give a full theoretical derivation of 
the coupled diffusive equation for spin dynamics in a fer- 
romagnetic metal film and describe the general form of 
the Rashba torque from weak to strong Rashba coupling 
limits. In Sec. |Hl we combine the Keldysh formalism 
and the gradient expansion technique to derive a diffu- 
sion equation for charge and nonequilibrium spin densi- 
ties. To demonstrate that the diffusion equation provides 
a coherent framework to describe the spin dynamics, we 
dedicate Sec. IHII to the spin diffusion in a ferromagnetic 
metal, which shows an excellent agreement to early result 
on the same system. In Sec lIVl we illustrate that the ab- 
sence of magnetism (in our diffusion equation) describes 
the well-known phenomenon of electrically induced spin 
polarization. The cases of a weak and a strong spin-orbit 
coupling are discussed in SeclVland Sec. [VH respectively, 
where we provide an analytical form of the Rashba torque 
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FIG. 1. (color online). A ferromagnetic thin film (such as Co) 
is sandwiched between a heavy metal layer (such as Pt) and 
an oxide (such as AlOx). A charge current is flowing in the 
ferromagnetic film along the x direction. 



where all Green's functions are full functions with inter- 
actions taken care of by the self-energies E R > A > K . The 
retarded (advanced) Green's function in momentum and 
energy space is 



G R{A \k,e) 



1 



e-e fe -&■ b(k) - £«( A )(fc,e) 



(3) 



where e k = k 2 /(2m) is the single particle energy. We 
have introduced a fe-dependent total effective held b(k) — 
A xc m/2 + a(k x z) of the magnitude bk — \A xc m/2 + 
a(k x z)\ and the direction b = b(k)/bf.. 

Neglecting localization effects and electron-electron in- 
teractions, we assume a short-range (S-function type im- 
purity scattering potential. At a low impurity concen- 
tration and a weak coupling to electrons, a second-order 
Born approximation is justified^ i.e., the self-energy due 
to impurity scattering is 2 ^ 

± R > A > K (r, r') = S ±lllG R < A ' K (r, r), (4) 



in an infinite medium. Section [VIII is dedicated to, for a 
wide range of relative strength between Rashba coupling 
and the exchange splitting, numerical solutions of the dif- 
fusion equation. In Sec. I Villi we discuss the implication 
of the Rashba torque on magnetization dynamics. 



II. DIFFUSION EQUATIONS 

As depicted in FigJTJ the system of interest is defined as 
a quasi-two-dimensional ferromagnetic metal layer rolled 
out in the x — y plane. Two asymmetric interfaces pro- 
vide a confinement in the z direction, along which the po- 
tential gradient generates a Rashba spin-orbit coupling. 
Therefore a single particle Hamiltonian for an electron of 
momentum fe is (K = 1 is assumed throughout) 

fc 2 - 1 

H = h an ■ (fc x z) + -A xc a m + H l (1) 

2m 2 

where <r is the Pauli matrix, m the effective mass, 
and m the magnetization direction. The ferromag- 
netic exchange splitting is given by A xc and a repre- 
sents the Rashba constant (parameter). The Hamilto- 
nian H % — J2j=i — r j) sums the contribution of the 
non-magnetic impurity scattering potential V(r) local- 
ized at rj. 

To derive a diffusion equation for the nonequilib- 
rium charge and spin densities, we employ the Keldysh 
formalism^ Using the Dyson equation, in a 2x2 spin 
space, we obtain a kinetic equation that assembles the 
retarded (advanced) Green's function G R (G A ), the 
Keldysh component of the Green function G K , and the 
self-energy T, K , i.e., 



where the momentum relaxation time reads 



[G 



R]-1qK 



G K [G A }~ 1 = ± K G A -G R ± K , (2) 



- « 2?r 

T 



d 2 k' 



\V(k-k')\ 2 5{e k -e k ,). 



(5) 



V(k) is the Fourier transform of the scattering potential 
and the magnitude of k and fc' is evaluated at Fermi 
vector hp. 

The quasiclassical distribution function g = gk,e(T, R), 
defined as the Wigner transform of the Keldysh function 
G K (r, t; r', t'), is obtained by integrating out the relative 
spatial-temporal coordinates while retaining the center- 
of-mass ones R = (r + r')/2 and T = (t + t')/2. As the 
spatial profile of the quasiclassical distribution function 
is smooth at the scale of Fermi wavelength, we may ap- 
ply the gradient expansion technique on Eq. ((5]) r§ which 
gives us a transport equation associated with macro- 
scopic quantities. Under the gradient expansion, the left 
hand side of Eq.([2]) becomes 

[G R ]- 1 G K ~G K [G A Y 1 

w [g,<r-b{k)]+ l -g + i^L 



i ( k 



a{z x &),V R g\ , (6) 



2 m 



where {•, •} is the anticommutator. The relaxation time 
approximation indulges the right hand side of Eq.Q as 

t K G A ~G R t K 

1 r /5(e, T, R)G A (k, e) - G R {k, e)p(e, T, R)] (7) 

where we have introduced the density matrix by integrat- 
ing out fc' in g, i.e., 



p(E,T,R) = 



1 



2vr^ n 



d 2 k' 



9k',e(T, R). 



(8) 



3 



For the convenience of discussion, the time variable is 
changed from T to t. At this stage, we have a kinetic 
equation depending on p and g 



i[& ■ b(k),g] + -g+^- + \\-+a{zx &),V R g 
t at 2 1m 



G R (k,e)p(e)-p(e)G A (k,e) 



(9) 



A Fourier transformation on temporal variable to the fre- 
quency domain u> leads to 



£lg - b k [Uk,g] = iK, 



(10) 



where = to + i/r and the operator Uk = cr ■ b satisfies 
UkUk = 1. The right hand side of Eq. fTO)) is partitioned 
according to 



K = 



1 f fc 

-- < — +a(z x <t), V R g 

2 m 



G R (k,e)p(e)~p(e)G A (k,e) 



(11) 



The equilibrium part is denoted by and the gradient 
term is treated as perturbation. Functions g and p 
are both in frequency domain. 



We solve Eq. (| 10[) formally to find a solution to g 



(2b 2 - n 2 )K + 2b 2 U k KU k - nb k [U k ,K] 



n(4bl 



n 2 ) 



J?[K]. 
(12) 



An iteration procedure to solve Eq. (TT2l has been outlined 
in Ref. |25| . We adopt the procedures here: according to 
the partition scheme of K, we use to obtain the 

zeroth order approximation given by g(°) = Jz? [K^(p)] 
which replaces g in to generate a correction due 

to the gradient term, i.e., (g^); we further in- 
sert K^(g^) back to Eq. (fT2"]) to obtain a correction 
Jt? [K^ 1 * (g( ')], then we have the first order approxima- 
tion to the quasiclassical distribution function, 



^=5 (0) +^ (1) (5 (0) )]. 



(13) 



The above procedure is repeated to desired orders using 

.^f^W^C"- 1 ))]. (14) 



3<»> =$c»-i> 



In this paper, the second order approximation is suffi- 
cient. The full expression of the second order approxima- 
tion for g is tedious thus to be excluded in the following. 
The diffusion equation is derived by an angle averaging 
in momentum space, which allows all terms that are of 
odd order in fej (i = x, y) to vanish while the combina- 
tions such as fcjfej contribute to the averaging by a factor 
kpSij 

Furthermore, a Fourier transform from frequency do- 
main back to the real time brings a diffusion-like equation 
for the density matrix, 



d 11 

j-p(t) =DV 2 p p + (z x cr) • p(z X cr) + iC [z X <r, V p] - B {£ x cr, V/3} 

+ T [(m x V) z p — a- ■ mV/5 • (z X cr) — (z x <x) ■ V per ■ m] 

+ (o - ' m P cr -m-p)- iA xc [& ■ m,p] - 2R{& ■ m, (m x V) z p} , (15) 



where all quantities are evaluated at Fermi energy cf- 
In a two-dimensional system, the diffusion constant D = 
tv f /2 is given in terms of Fermi velocity vf and mo- 
mentum relaxation time r. The renormalized exchange 



splitting reads A a 



(A, c /2)/(4£ 2 + 1) where £ 2 



(A 2 c /4 + a 2 k 2 F )T 2 . The other parameters are given by 



C 

R 
1 



(4e + 



I) 2 

^2 



aA xc v F k F T 2 
2(4£ 2 + l) 2 ! 



2(4£ 2 + 1) 



B = 



2a 2 k%T 



4£ 2 



1 



1 

T X C 



4£ 2 
A 2 



4£ 2 



1 



t xc is the relaxation time due to the D'yakonov-Perel 
mechanism^ Equation (|15|) is valid in the dirty limit 
£ <C 1, which enables the approximation 1 + 4£ 2 « 1. 
Charge density n and the nonequilibrium spin density S 
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are introduced by the vector decomposition of the den- 
sity matrix p = n/2 + S ■ <r. In real experiments ; 13 ' 19 ! 21 
spin transport in a ferromagnetic film suffers from ran- 
dom magnetic scatterers, for which we introduce, phe- 
nomenologically, an isotropic spin- flip relaxation S/r s f. 

Eventually, we obtain a set of diffusion equations for 
the charge and spin densities, i.e., 

dn 

— = DV 2 n + BV Z ■ S + TV Z • ran + RV Z ■ m(S ■ m), 

at 

(16) 



m x (S x m) 
T 

-L xr. 



and 

^. =DV 2 S _^l_S ± _ AxcSxm 

Ot T|| Tj_ 

+ BV z n + 2CV Z x S + 2R{m ■ V z n)m 

+ T[m x (V 2 x S) + V z x (m x S)} , (17) 

where V z = z x V. The spin density Sm = S x x + S y y is 
relaxed at a rate 1/tn = 1/t xc + l/r s f while S± = S z z 
has a rate l/rj_ = 2/t xc + l/r s /. 

For a broad range of the relative strength between 
the spin-orbit coupling and the exchange splitting, 
i.e., akp/A xc , Eq. (111)1) and Eq. flTTl) describe the spin dy- 
namics in a ferromagnetic film. When the magnetism 
vanishes {A xc — 0), the B term provides a source that 
generates spin density electricallyi 2 ^ On the other hand, 
when the Rashba spin-orbit coupling is absent (a = 0), 
the first two lines in Eq. (fl7l) describe a diffusive motion 
of spin density in a ferromagnetic metal, which, to be 
shown in the next section, agrees excellently with early 
results^ The C term describes the coherent precession of 
the spin density around the effective Rashba field. The 
precession of the spin density (induced by the Rashba 
field) around the exchange field is described by the T 
term, is thus at a higher order (compared to C) in the 
dirty limit for T = A xc tC/2. The R term contributes to 
the magnetization renormalization. 



III. SPIN DIFFUSION IN A FERROMAGNET 

Spin diffusion in a ferromagnet has been discussed ac- 
tively in the field of spintronics^ 7 . - - In this section we 
show explicitly that, by suppressing the Rashba spin- 
orbit coupling, Eq. fll~7|) is able to describe spin diffusion 
in a ferromagnetic metal. 

A vanishing Rashba spin-orbit coupling means a = 0, 
Eq.(17l) reduces to 



dS DV 2 S+ mX ^ ^ m x (S x m) 

8t TA T s f T xc 



(18) 



where ta = 1/A XC sets the time scale for the coherent 
precession of the spin density around the magnetization. 
This equation only differs from the result of Ref . [28| by a 
dephasing term of the transverse component of the spin 
density that is of the time scale T xc . 



In a ferromagnetic metal, we may divide the spin den- 
sity into a longitudinal component that follows the mag- 
netization direction adiabatically, and a deviation that is 
perpendicular to the magnetization, i.e., S = s^m + SS 
where sq is the local equilibrium spin density. Such a par- 
tition, after restoring the electric field by V — >■ V + eEd e , 
gives rise to 



= s DV 2 m 
SS 



DV 2 5S 

SqITI 



DeP F M F E ■ Vm 



- — + A xc m x SS, 

J- T.r. 



(19) 



where the magnetic order parameter is allowed to be spa- 
tial dependent m = m(r,t). The energy derivative is 
treated as d e S ~ PfMft^i given Pp the spin polarization 
and Nf the density of state, both are at Fermi energy 
e F - 

In a smooth magnetic texture, the characteristic length 
scale of the magnetic profile is much larger than the 
length scale for electron transport, we discard the contri- 
bution DW 2 5S^- The diffusion of the equilibrium spin 
density follows soDW 2 m w som/T s f. In this paper, 
we retain only terms that are at first order in temporal 
derivative, which simplifies Eq. (|19p to 

SS mxSS 



<T = DeP F N F E ■ Vm - 

TA T A 

The last equation can be solved exactly 



*£. (20, 



SS 



TA 



l + <r 2 



Pf 



m x (j e ■ V)m 



Pf,. 

S (je 

e 



-SQm x 



dm 



dm 



(21) 



where <; = ta(1/t s / + 1/T XC ) and the electric current 
j e = e 2 nTE/m is given in terms of electron density n. 
Apart from the transverse dephasing time implemented 
in parameter the nonequilibrium s pin density Eq.(|21j) 
agrees excellently with Eq.(8) in Ref. |27| . 

Given the knowledge of the nonequilibrium spin den- 
sity, the spin torque, defined as 



T = mx8S + — 

ta T xt 



SS, 



(22) 



is given by 
„ 1 



1 + ^ 2 



dm dm 

dt dt 



+ V~(3e 



V)m 



ft — m X (Je 

e 



V)m 



(23) 



where r\ = 1 + sta/T xc and f3 = T/\/r s f. Assuming a long 
dephasing time of the transverse component (i.e., T xc ~^ 
oo), then r\ k 1 and Eq. (1231) reproduces the Eq.(9) in 
Ref. [13] ■ On the other hand, a short spin dephasing time 
(i.e., T xc -> 0) yields 77/(1 + s 2 ) -> 1 and /3/(l + c: 2 ) -> 
which results in a pure adiabatic torque, i.e., the torque 
reduces to the first and third terms in Eq. 
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IV. ELECTRICALLY GENERATED SPIN 
DENSITY 

The effect of an electrically generated nonequilibrium 
spin density due to spin-orbit coupling 2 - can be ex- 
tracted from Eq. (fl7|) by setting exchange interaction zero 
(i.e., A xc — 0). Retaining D'yakonov-Perel as the only 
spin relaxation mechanism and letting r s / = oo, Eq. (|17p 
reads 



DV Z S 



S + S z z 



2CV z xS + 5V 2 n = (24) 



which reduces to describe the well-known spin-Hall 
effect i 25 i 31 i 32 Besides the spin relaxation [the second term 
in Eq. ([24]) ]. the spin dynamics is controlled by two 
competing effects: the spin precession around Rashba 
field (third term) and the electrical spin generation first 
pointed out by Edelstein. 2 In the case of an infinite 
medium along transport direction, i.e., the x direction, 
Eq. ({2"4")) leads to a solution 



S =eET xc B — y 



eEC, r 
ttvf ' 



-y- 



(25) 



where only the linear term in electric field has been re- 
tained. On the right hand side, we have used the charge 
density in a 2D system n — k 2 F /(2it) and introduced the 
parameter £ = ak F r as used in Ref . [25jj . 

It is worth pointing out that in the case of a weak 
spin-orbit coupling, only the spin precession term sur- 
vives, whereas in the case of strong coupling, the elec- 
trical spin generation dominates. As will be exposed in 
the remaining of this work, this distinction defines two 
regimes where the spin torque symmetry and efficiency 
are different. 



V. WEAK SPIN-ORBIT COUPLING 

A weak Rashba spin-orbit coupling implies a low 
D'yakonov-Perel relaxation rate 1/t xc oc a 2 , such that 
t~xc 3> T s f , ta , which allows spin relaxation to be dom- 
inated by random magnetic impurities. In this regime, 
when comparing the contribution due to C and T term 
to that from B and R, the later are at a higher order in 
a, thus to be disregarded. As mentioned in the previous 
section, this is a regime when spin precession about the 
Rashba field dominates the electrical spin generation. 

We consider a stationary state meaning dS/dt = 0. An 
electric field is applied along the x direction E — Ex. In 
an infinite medium^ all the spatial derivatives vanishes 
(V — > 0) and the dynamic equation reads 



2eECy x d e S + eET [y x (m x d e S) + mx(yx d e S)\ 

m x S m x (S x m) S 
= + ^ L + — 

TA J-xc T s f 



(26) 



In addition to spin density induced by the exchange split- 
ting, a weak spin-orbit interaction contributes to a devi- 
ation (in spin density) that can be considered as pertur- 
bation. Therefore, we may divide S = S± + Suva into a 
longitudinal and a transverse component. Equation (|26|) 
reduces to 



— m x S i - J— Si 



1 



Sum 



TA T± T s f 

= -eEP F Af F [2Cy x m + Tm x (y x m)], (27) 

where 1 /Tj_ = 1 /T xc + 1 /r s / and we have again employed 
the approximation d e S w PfA'f'TI and replaced the en- 
ergy derivative of the charge density by the density of 
states at Fermi energy, i.e., d t n w ti/ef — A/jr- We solve 
Eq. (|27p to obtain a nonequilibrium spin density 



S± = 



TA 



eEP F M F [(2C + ^r)m x (y x m) 



1 + ^ 2 
(T-2cC)(» x m) 



(28) 



and S\\ = 0. In Eq. ([28)) . the second component, oriented 
along the direction yxm, is actually perpendicular to the 
plane spanned by the magnetization m and the effective 
Rashba field that is along y. This field, as to be shown 
below, contributes to a torque that fulfils the symmetry 
described in recent experiments. 19 

The definition Eq. ([22]) leads to a general expression for 
the Rashba torque 



T =Tx_y x m + T»m x (y x m), 



(29) 



which consists of both out-of-plane and in-plane compo- 
nents with magnitudes determined by 



Ti = 



Til = 



eEP F M F 

1 + ^ 2 
eEPpAfp 

1 + ^ 2 



(2 V C + /3T), 
( V r-2(3C). 



(30) 
(31) 



Note that the sign of the in-plane torque, Eq. (|3Tj). can 
change depending on the competition between spin re- 
laxation and precession. 

To compare directly with the results in Ref. [l(J , we al- 
low the spin relaxation time r s f — > oo, then (3 w 0. We 
also consider the transverse dephasing time to be negligi- 
ble (T xc — > 0). Under these assumptions, 77/(1 + <r 2 ) ~ 1 
and we have T ± k 2eEP F M F C and Tu « eEP F Af F T. 
In the dirty limit, r -C C due to A xc t -C 1. By mak- 
ing use of the relation for the polarization P F = A xc /e F 
and the Drude relation j e = e 2 nrE /m, we obtain the 
out-of-plane torque 



amA xc . „ 

T = 2 j e y x m, 

ee F 



(32) 



which agrees excellently with the spin torque in an in- 
finite system in the corresponding limit as derived in 
Ref. 101. 
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VI. STRONG SPIN-ORBIT COUPLING 

The opposite limit to Sec|V]is a strong spin-orbit cou- 
pling. In this case, we consider that akp 3> A xc and 
the D'yakonov-Perel relaxation mechanism is dominat- 
ing, i.e., 1/t xc ^> l/r s f, since l/r xc oc a 2 . The spin pre- 
cession about Rashba field becomes negligible compared 
to the electrical spin generation (C << B) and therefore, 
it is not physical to simply assume that the direction of 
spin density is dominantly aligned along the magnetiza- 
tion direction, as what is treated in the case of a weak 
spin-orbit coupling. For the nonequilibrium spin density, 
a self-consistent solution from Eq. fTT)) is more justified. 

As in Sec|Vl we consider an infinite system where an 
electric field E is applied along the x direction, while 
the magnetization direction is left arbitrary. We approx- 
imate the energy derivative by d e ~ 1/cf- The above 
assumptions simplify Eq. flTH) to 



S x m m x (S x m) 2eEC 



1~A 



T, 



yxs 



{S + S z z) _ eEn 



By, (33) 



where a strong spin-orbit coupling renders Y and R terms 
negligible. By considering T xc ^> ta,t xc , Eq. (l33|) re- 
duces to 



S x m (S + S z z) 



2eEC „ „ e£n „ , niS 
-yxS=—By, (34) 



£f 



which is a set of linear equations for the nonequilibrium 
spin density. We are interested in the linear response 
regime, which implies that at the distance as defined by 
the Fermi wave length 1/kp, we have eE/kp -C akp. 
Therefore, up to the first order in exchange splitting, we 
extract the spin density from the above equation 



S = eEr xc B- 



y-XV xm 



X 



(35) 



where \ = t xc /ta and we have used the identity yxm = 
m z x — m x z. It is interesting to point out that, when 
taking the limit A xc — > hence x — > Oj Eq. (|3"S"|) precisely 
reduces to the spin density that is generated entirely due 
to spin-orbit interaction, i.e., Eq. (|25[) . 

The spin density Eq.Q35p yields a spin torque 

_ amA xc , „ 
T = j e (yxm 



eep 



+X rn x (y x m ) — -^m x z x m ) . (36) 



This torque is different from the weak Rashba limit and 
has a strong implication in terms of magnetization dy- 
namics. When the magnetization is aligned perpendicu- 
lar to the film plane m = z, 



1!) 



'-(m=z) 



amA xc 
eep 



j e (yxz + X zx(yxz)), (37) 



where the in-plane torque (second term) can be under- 
stood as either being generated from an effective field 
along the x direction or from a spin current that is po- 
larized along the y direction. In the strong spin-orbit 
limit, simply increasing the magnitude of Rashba cou- 
pling is not the most optimal way to enhance the in- 
plane torque efficiency since the D'yakonov-Perel spin 
relaxation rate increases as a 2 . The implication of this 
torque on the current-driven magnetization dynamics is 
discussed in Sec lVIIll 

In this paper, our numerical results as presented in 
the following section is based on the setup where m = x. 
The torque then reads 

T(m=a) = — 3e [y x x + -x x (y x x)j . (38) 

For the same parameters, an in-plane magnetization pro- 
duces an in-plane torque with a magnitude that is only 
half of the one with an out-of-plane magnetization, see 
Eq.fl37|). 



VII. NUMERICAL RESULTS 

We numerically solve the Eqs. p^j) and (IT71) to demon- 
strate that they provide a coherent framework to describe 
the spin dynamics in the diffusive regime for a wide range 
of parameters. Here, we consider an in-plane magnetiza- 
tion that lies along the x direction and we refer the read- 
ers to Ref. 12] for another case where the magnetization 
is perpendicular to the thin film plane. Using the finite 
element method, for a two-dimensional electron system 
we adopt the following boundary conditions: (i) vanish- 
ing spin accumulation at the edges along the transverse 
direction i.e., S(y — 0,L) = 0; (ii) an electric field is 
implemented along the x direction therefore we set the 
charge densities at two ends of the propagation direction 
to be constant np = np = np. The first boundary con- 
dition implies a strong spin-flip scattering at the edges, 
which is consistent with the experimental observations 
in spin-Hall effect £ The second boundary condition sets 
the charge density at the Fermi level. Equivalently, one 
can apply a voltage drop along the transport direction 
instead of an explicit inclusion of an electric field. 

The numerical results of nonequilibrium spin densities 
are summarized in Figj2j When viewing from the top 
panels (a,b) to the lower ones (c,d), for a fixed value 
of the exchange splitting, the system transits from the 
weak spin-orbit coupling regime to the strong spin-orbit 
coupling regime. As an illustration of this transition, 
the S z component of the spin accumulation evolves from 
a symmetric spatial distribution in the weak spin-orbit 
coupling regime [a = 0.0005 eV nm in Fig. 1(a)] to an 
antisymmetric spatial distribution along y in the strong 
coupling regime [a = 0.05 eV nm in Fig. 1(c)]. Note that 
throughout this transition, the in-plane spin density S y 
is robustly constant in the bulk of the wire. 



7 




-25 -15 -5 5 15 25 -25 -15 -5 5 15 25 



>' (nm) y (nm) 

FIG. 2. (color online). Spatial profile of the nonequilibrium 
spin density S z (a),(c) and S y (b),(d) for various values of the 
Rashba constant. The width of the wire is L = 50 nm. The 
magnetization direction is along the x axis. Other parame- 
ters are: momentum relaxation time r = 10~ 15 s, exchange 
splitting A xc — 0.01 eV, spin relaxation time r 3 / = 10~ 12 s, 
and the Fermi vector k,F = 4.3 nm~ . 



x 10 x 10 




a (eV nm) a (eV nm) 

FIG. 3. (color online). The magnitude of the out-of-plane 
torque T± (a) and in-plane torque Tji (b) as a function of 
Rashba constant for various exchange splitting. Other pa- 
rameters are the same as in Fig[2] 

a decrease in the torque magnitude. In fact, the transi- 
tion suggests that the optimal magnitude of the in-plane 
torque is achieved when the exchange energy is about the 
same order of magnitude of the Rashba splitting akp- 



This symmetry transition and the emergence of peaks 
close to the boundaries are a manifestation of the com- 
petition of the Rashba and exchange fields. In the weak 
coupling regime, the total field is dominated by the ex- 
change field pointing at the x direction, around which the 
spin accumulation profile that is spatially symmetric. As 
the spin-orbit coupling increases, the total field is tilted 
towards the y axis, then the spin projections along +y 
and — y are no longer symmetric, as indicted by curves 
with intermediate a values in Figj^a)(b). In the strong 
coupling regime, when Rashba coupling is larger than the 
exchange field, the antisymmetric profile in S z and the 
symmetric one of S y follow naturally from the spin-Hall 
effect induced by the Rashba spin-orbit interaction. 

The out-of-plane and in-plane torques are plotted in 
Fig 13] with respect to the Rashba constant a for various 
exchange splitting. The transition regions are of par- 
ticular interest. During the transition from the weak to 
strong coupling, the magnitude of the out-of-plane torque 
T± (see Fig|3][a)) first reaches a plateau, then rises again 
as the a increases. In the large a limit, though the mag- 
nitude of the torque increases with a, the torque effi- 
ciency defined as dT±/da is actually smaller than it is 
in the weak coupling. This picture is consistent with a 
semi-classical Boltzmann equation description employed 
in Ref. [l(|. This behavior is due to the different pro- 
cesses generating the Rashba torque in both regimes. As 
discussed in Sec. V and VI, in the weak coupling regime, 
the torque is dominated by the spin precession around the 
Rashba field; whereas in the strong coupling, the electri- 
cal generation of the spin density dominates. These two 
distinct processes have different efficiencies. 

The in-plane torque Tj| behaves differently. In the 
strong coupling limit, Ty is proportional to 1/a due to the 
large D'yakonov-Perel spin relaxation rate that is of the 
order a 2 . Therefore a stronger spin-orbit coupling means 



VIII. DISCUSSION 

Current-induced magnetization dynamics in a sin- 
gle ferromagnetic layer has been observed in vari- 
ous structures that involve interfaces between transi- 
tion metal ferromagnets, heavy metals and/or metal- 
oxide insulators. Existing experimental systems are 
Pt/Co/AlO^^i&ZI Ta/CoFeB/MgO^ Pt/NiFe and 
Pt/Co bilayers.— Besides the structural complexity in 
such systems, an unclear picture of spin-orbit coupling 
in the bulk and at the interfaces places a challenge to 
understand the nature of the torque. 



A. Validity of Rashba model 

The celebrated Rashba-type effective interfacial spin- 
orbit Hamiltonian was pioneered by E. I. Rashba to 
model the influence of asymmetric interfaces in semicon- 
ducting two-dimensional electron gas£ a sharp potential 
drop, emerging at the interface (say, in the x — y plane) 
between two materials, gives rise to a potential gradient 
W that is normal to the interface, i.e., VV « £,{r)z. In 
case a rotational symmetry exists in-plane, a spherical 
Fermi surface assumption allows the spin-orbit interac- 
tion Hamiltonian to have the form Hr = ever ■ (p x z), 
where a ss (£)/4m 2 c 2 . As a matter of fact, in semicon- 
ducting interfaces where the transport is described by a 
limited number of bands around a high symmetry point, 
the Rashba form can be recovered by k ■ p theory^ 

As far as metallic interfaces are concerned, a spin-orbit 
splitting of the Rashba-type in the conduction band has 
been observed at Au surfaces Gd/GdO interfaces^ Bi 
surfaces and compounds, 36 and metallic quantum wells^l 
The presence of a Rashba interaction at oxide hetero- 
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interfaces 3 -^ has also been reported recently. It is worth 
noticing that the symmetry breaking-induced spin split- 
ting of the conduction band seems rather general and 
might not be restricted to heavy metal interfaces^ 

In transition metals, however, the free electron approx- 
imation fails to accurately characterize the band struc- 
ture due to both a large number of band crossing at 
the Fermi energy and a strong hybridization among s, p 
and d orbitals. Density functional theory (DFT) is suc- 
cessful in investigating spin-orbit interaction at metallic 
surfaces! 39 ' 40 For example, in Refs.[40|. the authors ob- 
serve a band splitting that possesses similar properties as 
Rashba spin-orbit interaction and decays exponentially 
away from the surfaced Alternatively, the spin-orbit in- 
teraction at metallic surfaces has been addressed using 
tight-binding models for the p orbitalsi 41 i 42 At such sharp 
interfaces, the magnitude of the orbital angular momen- 
tum (OAM) is considered to play a dominant role at the 
onset of a Rashba-type spin splitting. 

This finding is consistent with the long standing 
work on interfacial magnetic anisotropy at a ferro- 
magnet/heavy metal^ 3 - and ferromagnetic/metal-oxide 
interfaces^ where a perpendicular magnetic anisotropy 
arises from the orbital overlap between the 3d states of 
the ferromagnets and the spin-orbit coupled states of the 
normal metal. The observation of perpendicular mag- 
netic anisotropy at Co/metal-oxide interfaces tends to 
support the major role of large interfacial OAM in the 
onset of interfacial spin-orbit effects i 42 ' 44 

The presence of interfacial Rashba spin-orbit coupling 
has also been shown to explain tunneling anisotropic 
magnetoresistance 4 ^ and interfacial perpendicular mag- 
netic anisotropy— Assuming a simplistic free electron 
dispersion, these models properly capture the symmetries 
of the observed effects, demonstrating the seminal role of 
interfacial Rashba spin-orbit coupling in these systems. 

All the previous theoretical and experimental studies 
strongly suggest that the interfacial spin splitting exists 
in the presence of a large OAM and potential gradient. 
However, a microscopic description of realistic interfaces 
is absent. Although the Rashba spin-orbit interaction is 
a convenient Hamiltonian to extract qualitative behav- 
iors, its applicability to realistic metallic interfaces with 
complex band structures remains to be tested. Mean- 
while, a recently proposed phenomenon called spin-orbit 
sound may provide an alternative experimental tool to 
determine the strength of the coupling.— 



(jf(SH) ^ ue ^ Q g p m _jjall effect) that is transverse to the 
magnetization direction* 2 ^ In the bulk, the spin current 
can be estimated using the ratio between spin- Hall {<J SH ) 
and longitudinal {a xx ) conductivities (the so-called spin- 
Hall angle), i.e., 



3 



(SH) _ 



-Je- 



(39) 



A perturbation calculation using the second-order Born 
approximation gives rise to a spin current thus the torque 
with a magnitude given by, 



T (SH) = 



2er° : 



(40) 



where, in general, 7 > 1 is a dimensionless parameter 
taking into account both side-jump (7 = 1) and skew- 
scattering (7 > 1) contributions to the spin-Hall effect. 49 
r] so is the spin-orbit coupling parameter and r t ° r is the 
transport relaxation time due to bulk impurities, the 
same definitions as in Ref.[49| except here the definition 
of spin current differs by a unit l/(2e). Meanwhile, the 
magnitude of the in-plane torque Eq. pip can be simpli- 
fied to, since A xc t <C 1 



Tn 



4am 



-Je 



(41) 



The spin-orbit coupling parameter a in our definition in 
Eq.fll]) has the unit of energy. Equations (|4Tj) and (|40|) 
actually show that the in-plane Rashba torque and the 
spin-Hall torque have a very similar parameter depen- 
dence. 

Using a diffusive description of the bilayer system con- 
sisting of fcrromagnet/heavy metal, it was shown that 
both SHE torque and Rashba torque adopt a similar 
form, 



T = T»m x (y x m) + T±_y x m. 



(42) 



The similarity in the geometrical form of the two torques 
imply that, in principle, they are able to induce the same 
type of magnetic excitation. 23 



C. Magnetization Dynamics 



B. Comparison between SHE torque and Rashba 
torque 

At this stage, it is interesting to compare the parameter 
dependence of the in-plane torque Tji (in Eq. (|3"Tj) ) and the 
torque generated by a spin-Hall effect^ in the bulk of a 
heavy metal material such as Pt. In the later case, the 
torque T^ SH ' exerted on the normal metal/ferromagnet 
interface is obtained by projecting out the spin current 



In Pt/Co/AlOz trilayers, Miron et al have observed 
a current-driven domain wall nucleation,— an enhanced 
current-driven domain wall velocity^ 1 - and a magneti- 
zation switching^ The symmetry of the spin torque 
needed to explain the experimental findings agree well 
with Rashba torque proposed in Ref. [13 and [l2[ • In a 
similar structure, Pi et al^ and Suzuki et al^ also ob- 
served an effective fieldlike torque to be interpreted by 
the Rashba torque. 
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1. Magnetization switching 



integrating over the magnetic volume 



Our previous discussions suggest that both Rashba 
torque and SHE torque possess a general forrn^ 



T = T±y x m + T\\m x (y x m) 



(43) 



that consists of a fieldlike torque (T±) and an (anti- 
)damping term (Tn). As a consequence, both Rashba 
torque and SHE torque have the appropriate symmetry 
to excite the magnetization of a single ferromagnet and 
induce switching, as observed by Miron et a& and Liu 
et al£2- In the case of a large Rashba spin-orbit coupling, 
the torque acquires an additional component that acts 
like an effective magnetic field along the z direction, van- 
ishing as the magnetization component m x is zero (see 
Section IVI|) . which provides an additional torque that 
helps destabilize the magnetization. 



2. Current-driven domain wall motion 

The influence of Rashba/ SHE torque on a domain wall 
can be illustrated within the rigid Bloch wall approxima- 
tion. The perpendicularly magnetized Bloch wall is pa- 
rameterized by m — (cos <f> sin 6, sin cj> sin 9, cos 6) where 
(f) = cj)(t) and 6(x,t) = 2t&n- 1 [e (x - x ^/ A ], where x c 
refers to the center of the domain wall and A is defined 
as the domain wall width. To describe the dynamics of 
a Bloch wall, Landau-Lifshitz-Gilbert (LLG) equation 



dttn = — jm x H e ff + a G dtm x m + r (44) 

has to be augmented by the current induced torque t 

t = bj(u ■ V)m — (3bjm x (u ■ V)m 

+bj(r±y x m + T\\m x (j/x m) 

+ TzTTlxZ x m). 

(45) 

The torque r is written in the most general form, where 
the first two terms are the regular adiabatic and the so- 
called non-adiabatic torques; the next two terms (t\\ and 
t_l) emerge from the presence of Rashba and/or spin 
Hall effect and the last term r z appears only in large 
Rashba limit (see Sec. I VI I) . The magnitude of the adia- 
batic torque is bj = [iBPje/e- The effective field is given 
by 

2A 

H eff = ^-V 2 m + H K m x x + Hj_m z z. (46) 

Parameter 7 in LLG is the gyromagnetic ratio, a G is the 
Gilbert damping, A is the exchange constant, M s is the 
saturation magnetization, Hk is the in-plane magnetic 
anisotropy and H±_ is the combination of an out-of-plane 
anisotropy and a demagnetizing field. The magnetization 
dynamics can be obtained readily from Eqs. P4]) - (|4^1) by 



d t 4> + a G 



d t x c 

~7T 

d t Xc 



Ait 
— (-II 

-7— 



2 

sin 2(f) 



) cos <fi — f3 
1 



bj 
A 



(47) 
bj 



An 

— — Tj_ cos „ , 
2 J A 



(48) 



We observe that the in-plane torque T|| distorts the do- 
main wall texture, while the perpendicular torque t± 
drives the domain wall motion. The additional torque 
t 2 , arising in the large Rashba limit, only contributes to 
the in-plane torque. Therefore, in the following, we will 
refer to the in-plane torque as T|* = T|| — r z /2. Below the 
Walker breakdown (dt<fi = 0), the velocity is given by 



P 



d t x c = - 
7— sin 20 



An 

a G 



OLG 



(49) 



An 



{a G T± 



Tit ) COS ( 



a G A 



(50) 



where the tilting angle <fi is given by the competition be- 
tween the magnetic anisotropy, the non-adiabatic torque, 
and the Rashba/SHE torque. In the case of weak Rashba 
(t z — 0), assuming tii — (3t± and omitting the correction 
to the spin precession, we recover the results of Ref. [HJ . 
When neglecting the in-plane torque and accounting for 
the perpendicular Rashba torque (T|* = 0), the Rashba 
torque only acts like an effective transverse field and en- 
hances the Walker breakdown limit^i [see Eq. (|50])]. 

Accounting for the in-plane component T|| arising ei- 
ther from corrections to Rashba torque or from the SHE, 
this torque appears to modify the domain wall veloc- 
ity. Therefore, depending on the strength and the sign 
of Rashba/SHE torque as well as on the resulting tilt- 
ing angle 0, it is possible to obtain a vanishing or even 
a reversed domain wall velocity, as has been shown nu- 
merically in Ref. [52[ and illustrated in Eq. (|49]) . A full 
scale numerical investigation is beyond the scope of this 
article, but it will help understand the profound effect of 
Rashba and SHE torque on the domain wall structures. 



IX. CONCLUSION 

Using Keldysh technique, in the presence of both mag- 
netism and a Rashba spin-orbit coupling, we derive a spin 
diffusion equation that provides a coherent description to 
the diffusive spin dynamics. In particular, we have de- 
rived a general expression for the Rashba torque in the 
bulk of a ferromagnetic metal layer, at both weak and 
strong Rashba limits. We find that the torque is in gen- 
eral composed of two components, a field-like torque and 
an (anti-)damping one. Being aware of the recent alter- 
native interpretation on the current-induced magnetiza- 
tion switching in a single ferromagnet, we have discussed 
the difference between the Rashba and the SHE torques. 
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While exploring the common features, we found that the 
magnetization dynamics driven by the Rashba torque 
presents several interesting similarities to that induced 
by SHE torque. Nevertheless, further investigation in- 
volving structural modification of the system is expected 
to provide a deeper knowledge on the nature of the inter- 
facial spin-orbit interaction as well as the current-induced 



magnetization switching in a single ferromagnet. 
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